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ABSTRACT 

This paper describee a method for correcting antenna 
temperatures observed with a single-reflector antenna for the 
effects of side lobes. The contribution from side lobes 
is a severe problem for earth-viewing microwave antennas » 
since the contrast in brightness between land and ocean 
is very high. The method developed here for correcting 
for the antenna side Icbes estimates each brightness temperature 
as a linear ocmibination of the antenna temperatures; 
therefore, it is suitable for conditions where the data 
rates are very high. 

Any linear combination of the observed antenna temp- 
eratures is equivalent to the antenna temperature that would 
have been observed by an aperture with some other 
antenna pattern, which we call the "effective" power pattern. 

It is shown that this effective power pattern can be chosen 

so that it has lower sldelobe levels than the original power 
pattern. 

Another problem is to match the areas observed at 
two frequencies (where the beam patterns will be dif- 
ferent) so that meaningful comparisons can be made. This 
can be done using a simple extension of the material reported 
here. 

A simple (Gaussian) beam pattern is chosen as an 
example; some representative results show that it is 
possible to suppress the side lobes at the expense of 
only minimal noise ^unpliflcation and broadening of the 
main beam. It is also possible to smooth the observations 
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so that the noise Is reduced, the side lobe level reduced, 
and the main beam broadened only by about lo%. 
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INTRODUCTION 

This paper describes a method for processing the data 
from a scanning microwave radiometer so as to modify the 
antenna pattern of the Instrument. This method can be em- 
ployed for correcting the data from mechanically-scanned micro- 
wave Instruments that view the earth from space. Two such 
Instruments— the Scanning Multlchaxmel Microwave Radiometer 
on NIMBUS-7 Md SEASAT-1 and the Large Antenna Multichannel 
Microwave Radiometer under study for the NOSS satellite — ^use 
a single antenna to make observations simultaneously at fre- 
quencies from 6.6 to 37 GHz. These antennas are conically 
scanned; l.e., the antenna points away from the nadir axis and 
Is rotated about the nadir axis. The result Is that the bore- 
sight traces a cone whose axis Is In the nadir direction. The 
angle of Incidence of the boreslght on the earth's surface Is 
constant. The two problems to be addressed In the data pro- 
cessing are: 1) to modify the antenna pattern at the lowest 

frequency to reduce the level of the side lobes / and 11} to 
smooth the higher-frequency observations so that they can be 
compared with the lowest-frequency observation. 

Many authors have studied the resolution limits Imposed 
by an aperture of finite size (e.g., Schell 1969, and Buck and 
Gustlnclc 1967) . Methods have been devised for correcting for 
side lobes In Interferometer measurements (Hogbom 1974) or 
for single-reflector emtennas (e.g., Chow and Pelletier 1974, 
Emerson et al. 1979, and Njoku 1980). All of these methods 
require extensive data processing; the data rates from satellite 
Instruments eure high enough to preclude the use of these methods. 


Also, attention seems to have been focused on trying to 
improve the resolution in the main beam (e. g. , Classen and 
Fung 1974 and Stogryn 1978) ; it is shown in a later section 
that this must be accompanied by amplification of the noise 
in the observations. The focus of this naoer is to consider 
ways to decrease the level of the side lobes of the antenna 
without amplifying the noise greatly. In general, little 
can be effectively done to improve the main- lobe resolution 
(because of noise amplification) , but the effects of the 
side lobes can be suppressed. 
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The method presented here is to find a set of coefficients 
so that the estimated brightness temperature at a point is 
given by 

where is the anteima temperature observed at the point (x^,y^). 

Let P(X/y) be the true power pattern of the ant.tnna , and 
Q(x,y) be a desired power pattern. The aim is to find a set 
of coefficients so that the estimated brightness temperature 

A 

Bj is the same as the antenna temperature that would have been 
observed with an antenna that had the power pattern Q(x,y). It 
will turn out that Q(x,y) cannot be chosen arbitrarily— e.g. , it 
is not possible to use a delta function — and that any increase 
in the spatial resolution will be accompanied by an increase in 
the noise in the estimated value B. In this paper, emphasis 
will be placed on finding patterns Q(x,y) that have low side- 
lobe levels. Suppression of the far sidelobe response in im- 
portant for earth-viewing satellites, since at microwave fre- 
quencies the contrast in brightness temperature between land and 
ocean is on the order of 100 K; observations of the ocean will 
be useless when a few percent of the power comes from land areas. 

Section I will discuss the mathematical methods employed to 
correct the antenna pattern using a matrix technique. A parti- 
cular antenna power pattern is chosen for study and described in 
Section II; some numerical results are given in Section III. 

In particular, estimates are found for the error in the estimated 
brightness temperature and the factor by which the instrument 
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noise is amplified in the data processing. The proolem of 
unwanted noise amplification is a common one in remote 
sensing problems — e.g.« see Foster (1961) and Fleming (1976) 
for a review of related problems. Methods that do not 
treat the noise amplification problem from the beginning- 
such as the one proposed by Njoku (1980) -—usually run into 
serious difficulties. 

In Section IV, the problem will be looked at by treating 
the antenna as a spatial filter (this is, in part, an exten- 
sion of the discussion given by Schell (1969) ; it also 
addresses a signific«mtly different problem) . This analysis 
is included to illustrate the properties and limitations 
of antenna pattern modification; it does not produce a work- 
able method for performing the required calculations. 

The problem of correcting the antenna patterns is thus 
addressed from two viewpoints. One is that of matrix oper- 
ations, where the elements of the matrices are integrals of 
the antenna pattern: this formulation is directly applicable 

to developing algorithms. The other viewpoint is that of 
analysis of the problem by treating the antenna as a spatial 
filter, and using Fourier transforms of the antenna pattern: 
this approach leads to an understanding of the problem, but 
is not suitable for performing calculations. 



I. MATHEMATICAL MODEL 


a) The matrix formulation 

In this section the matrix formulae used to perform the 
antenna pattern modifications are derived. Let P(x,y) be the true 
antenna pattern at the lowest frequency used for observations; the 
problem is to find a set of coefficients that transform P(x»y) into 
some more suitable effective pattern* We suppose that the antenna 
is pointing toward a plane far from the antenna; the plane is normal 
to the boresight of the antenna and the antenna is moved to various 
points (Xi,yi) . 


t^hen the antenna is pointing at , the antenna 

temperature is given by 
0 » 00 

»J J P^(x,y)B(x,y) dxdy (1) 

«00 *w00 

where B(x,y) is the brightness distribution and 

Pi(x,y) » P(x^-x,y^-y) . (2) 

The estimated brightness temperature B^ at the point (x^y^) 
will be given by 
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where » in fact, the coefficients may depend on the position 
(XjrYj)* 0(x,y) be the desired effective antenna pattern* 


Since 


OS OS 

Z)M^P^(x,y)B(x,y)dxdy 
•W i 


( 4 ) 


we want 


Q(x,y) -SM^P^Cx^y). 


( 5 ) 


The coefficients can be found from a least-squares fit of 
P(x,y) to Q(x,y) . 

The error in the fit will be 

• • f 2 

X* • 8* J J |a(x,y) -S-\Pi(x,y)j dxdy + 

<•00 -oo ^ I 


(6) 


where is the ncise in the radiometer, s^ the variance of 
B(x,y). 
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Note that consists of two terms r one is due to the 
difference between Q(x,y) and the fitted value. The error 
introduced by that term is orooortional to the variance in 
the brinhtness distribution s*. The other term is due to 
amplification of the noise in the measured antenna temperatures. 
The constMts are found from differentiating eguaticn € 
with resoect to each and equating the result to zero. 

The resu..:Jng effective anteuna pattern is given by 

'i'(x,y) ■ I M^P^(x,y). (7) 

The difference between equations 5 and 7 is that 0(x,y) is 
what is desired as the effective oower pattern; S!(x,y) 
is what can be achieved. 

Define 

OD 00 

Ri ■ II Pj^(x,y)Q(x,y)dxdy (8) 


and 




P^(x,y)Pj (x,y)dxdy. 


(9) 




Then the values of tfvat minimize x* are given by 


M ■ |p + (a^/s^)l| R 


( 10 ) 


where the elements of £ are . The matrix | is positive 
definite, so none of the eigenvalues of (P + (a^/s^)!) is 


less than or equal to zero. The matrix is non»singulir and 
its inverse exists. This solution is similar to one given 
by Foster (1961) . effect of using equation 10 is to 

balance the closeness of the fit of 4* to Q against letting 
the noise amplification become too large. 
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b) Pattern-matching at higher freauencles 

The method for matching the effective patterns at higher 
frequencies Is a simple extension of v^at has been done so far. 
Since the beam Is broadest at the lowest frequency » the 
effective pattern 4*^(x,y) must be found first for the 
lowest frequency. Then, letting 'Pj^(x,y) be the effective pattern 
at a higher frequency, the requirement that is 

eouivalent to finding a different set of coefficients 
that satisfy 

’?^(x,y) » H'j^(x,y) = I Pj^ (x-x^,y-y^) (12) 

where P^(x,y) Is the antenna pattern at the higher freauency. 

A solution to equation 12 always exists with J M/ * <1 

that corresponds to a smoothing of the antenna temperatures 
at the higher freouency. 
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II. THE ANTENNA PATTERN 

This section describes the antenna pattern used to illus- 
trate the method described in the previous section. This 
pattern was chosen so that the computations can be handled 
analytically. 

The antenna gain pattern is given by 

P(0) * 0.1475 + 0.0067 (1) 

where 9^ = x* + y*. P(0) is normalized so that its integral 

over all space is eaual to one. The approximate half-power 
beamwith (FWHM) is 2/Hn?. If the main beam is defined to 
be 2.5 times the FWHM, then the main-beam efficiency is 90%. 

This pattern has significant contributions from the side 
lobes (cf Figure 1-a) ; in the next section results are given 
for attempts to reduce those side lobes. 
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III. NUMERICAL STUDY 

Examples are given in this section of how the real power 
pattern of an antenna can be modified to suppress the side 
lobes. Some care must be taken In choosing the form for Q(x,y ) , 
the desired pattern. After some trial calculations. It became 
apparent that fitting to a boxcar function (e. g., eauatlon 
IV- 1) gives poor results; the main beam is fattened and the 
sidelobe levels are increased . A reasonable form for Q(x,y) 
is a Gaussian 


^ -0V2a* 


Q(x,y) * e 

2iTa^ 

where 6* * x* + y* . Then 


( 1 ) 


0.1475 0.0067 

R. = exp{-eJ/(2a*+2) } + ;; expi-6?/ (2a*+10) } (2' 


a^ + 1 


.2a*+ 1 


where 6? = x^ + y?. Then the coefficients are given by 


'■ij 


0.0494 


exp(-0?j/4) + 0.179 exp(-0?j/12) 


+ 0.0144 exp(-0|j/2O) 


(3) 


2 2 
where 6?^ = (x^ ” * ^^i ” * 
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If Nj^P^(x«y) Is going to represent an effective antenna 
pattern# it is necessary that 2^ ■ 1. However# there is a 

unique solution to equation X-9; in general# the coefficients 

will not quite be normalized. Therefore define c 

is a normalized set of coefficients. This is the solution to 
the equation 

(g, + a cR (5) 

= c Q(x,y); * o^/s^. The error 

be written 

+ c^K + n* 2 M?. (6) 

and substituting into equation 7# 

(7) 

Experience has shown that# if Q(x,y) is suitably chosen# c differs 
from unity by only a few percent. 

One important parameter in the modelling of the antenna pat- 
tern correction is the spacing of the observations of antenna 
ter.«perature Ax. Once the orbital altitude (and therefore the 
satellite velocity) and the scan rate have been chosen# the time 

spent observing a region is proportional to its area; therefore 

_2 

the noise in T^ is proportional to (Ax) . This effect has been 
included implicitly in the calculations reported here; the 


which arises if we let Q' (x#y) 
in the fit— eouation 1-6— can 

x’ - E - 2oEMj,R3 

Multiplying equation I>9 by M^ 
= S^(C*K - CM*R) 
where K *?/ Q(x#y)^dxdy. 
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quoted slgnal-to-noise ratio n” is divided by (Ax) i.e.» the 

signal-to-noise ratio is always given as the ratio of the variance 

in brightness to the noise for an integration time corresponding 

to Ax « 1. 

b) Results 

The problem now is to choose values of Ax, the interval 
at which the antenna temperatures are sampled , a‘ , the width 

of the function Q(x,y), and the signal->to-noise ratio rfS which 

is treated as a free peurameter. We have performed calculations 

with a variety of different values of these parameters; some 

combinations give good results, and others do not. It is not 

possible to quemtify each resulting antenna pattern with only 

a few numbers; some representative patterns are shown in Figure 1, 

and more are described in Tables 1 and 2. Since the purpose 

of modifying the antenna pattern is to reduce the level of the 

far sidelobes, these tables list the values of x for which 

the power is reduced by 10 db and by 20 db. Where there 

is a sidelobo level is higher then -20 db, the point «diere the 

main lobe is -20 db from the peak, and also the point where 

the side lobe is -20 db from the peak are both given. 

The levels of the side lobes are given in the last column. 

The first column of the tables lists the value of a^ , 

the second the signal-to-noise ratio used (rf^) , and the 

third the resulting noise amplification factor a^. 

Figure 1 shows the patterns when Ax « 1 (i. e., the 

antenna temperatures are sampled every half-beamwidth) , a^ » 1, 

2 3 

and Tf* = 10, 10 , 10 , and <». Figure 1-a shows the original 
antenna pattern and the pattern for a Gaussian of width 1. 

In some ways, the most interesting pattern is shown in 
Figure 1-b, where 10; the noise amplification factor 

is 0.079. That is, this pattern represents a considerable 
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smoothing of the raw data, but the effective antenna pattern 
Is not stronger at x > 4.2 than the original pattern. This 
pattern has. Instead, the property that the noise Is reduced, 
and only the main beam Is fattened. Figures 1-c to 1-e 
show what happens as the slgnal-to*nolse ratio Is Increased: 
the main beam gets narrower, the noise amplification Increases, 
and the sldelobes (which are alternately negative and positive) 
appear. Figure l~e shows that a Gaussian beam with width 1 
can be achieved with a noise amplification of only 1.58. 

Figure 1-f Is Included to show that It is not possible 
to narrow the antenna pattern beyond certain limits. 

This is an attempt to fit the antenna pattern to a Gaussian 
with a^ - <» . Even though the noise amplification 

factor Is 235, the fit Is not all that good: there are 
three sldelobes at about »20 db or higher. 

Figure l-g Is Included to show how the spacing of tl:e 

observed antenna temperatures Ax Influences the the width 

of the effective antenna pattern. The curves In Figures 

2 

1-c and 1-g both have S/N » 100 and a » 1; the sampling 
interval in Figure 1-g Is twice that of Figure 1-c. 

The former pattern is not an Improvement over the original 
pattern at all; the side lobes are higher and the main 
beam wider. In general, it seems from the computations we 
have done that the minimum sampling interval will be about 
one half of the half -power beamwidth (the maximum sampling 
interval size will be discussed in the next section) . 
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The quality of the corrected pattern can be measured 
by two additional numbers. The resolution is measured by 
the full width at half-maximum (FWHM) of the beam. A measure 
of how well the side lobes are suppressed is given by 
Let the brightness distribution be given by B ■ 1 for x < 
and B « 0 for x > x^. Then the antenna temperature is a function 
of Xj^; *1000 value of x^ where the antenna temperature 

< 10“^. For the original beam, the FWHM is 2.4, and x^qqq * 5.9. 
Tables 3 and 4 list the values for the corrected emtenna pat- 
terns. 

From Table 3, we can see that, for a* » 1, S/N » 1000, 
and Ax B the FWHM is 2.6 (10% larger than the original 
beam), but » 2. The noise amplification is 0.5. There- 

fore, for this combination of parameters, the distance at 
which the contribution from the side lobes » 10~^ is reduced 
from 5.9 to 2.0; the main beam is broadened only by 10%, 
and the noise is reduced by a factor of 2. This is certainly 
an improvement over the original beam. 
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IV. FOURIER TRANSFORMS AND 
THE RELATIONSHIP TO A LINEAR FILTER 

So far, the discussion of the antenna pattern correction 
problem has been in terms of integrals of the power pattern 
Md matrix operations. The following discussion will be in 
terms of the Fourier transform of the effective power pattern 
and the aperture illumination distribution. The discussion 
is limited in two respects that make it inapplicable to the 
specific problem of correction the data from a mechamically 
scanned radiometer looking at the earth, but the discussion 
will aid in the understanding of the processes. The restrictions 
are 1) the euitenna temperatures are observed at regularly- 
spaced intervals A, and 2) the power pattern does not change as 
function of position (for a conically-scaxming antenna, the 
power pattern rotates as the antenna scans) . 

In the previous discussion, the function Q(x,y) was intro- 
duced as the "desired" antenna pattern; there was, however, no 
indication as to what a desirable form might be for Q(x,y). 

In this section, the Fourier transform of Q(x,y) will be seen to be 
a smoothing function that broadens S'(x,y) (the effective power 
pattern) ; it turns out that this smoothing is what makes it possible 
to reduce the level of the sidelobes. 
a) Relation to aperture illumination distribution 

To start off, let P(x,y) be the power pattern, and 
P(x,y) * P(-x,-y). Let Q(x,y) and P(x,y) have Foiurier transforms 
q(s,t) and p(s,t); the Fourier transform of f(x,y) is p(s,t) . 

(Note that p(s,t) is the complex conjugate of p(s,t)). 
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The antenna temperatures are sampled at an interval d; 
define the function 


H(x) - 1 |x| < ij 

n(x) - 0 |x( > >s . 


( 1 ) 


The coefficients were defined in equation 1-8; 

is the convolution P (x^y) *P(x,y) evaluated at (x^-x^, 
yj^y^) (whore f*g denotes the convolution of f with g) . 

The Fourier transform of P*P is p(s«t)p(S/t) . The coefficients 
defined in equation 1-7 are values of Q(x,y) *P(x,y) and the 
Fourier transform is q(s,t)p(s,t) . If the expression for 
(equation 1-6) is differentiated with respect to and the 
result set equal to zero, the resulting equations can be written 


f _ 2ffi{(x.-x^)s + (Yi-y^ 

[p(s,t)p(s,t) J M^e ^ J ^ J 

2iri(x^s + y^t) 


-Y4)t} 


( 2 ) 


- p(s,t)q(s,t)e 


dsdt + » 0. 


Now introduce the function 

-2iri(x.s + y.t) 

m(s,t) « I MjC J J ; (3) 

this finite Fourier treuisform has the inverse 
7 7 2TTi(x,S + y.t) 

I ] n(s/A)H(t/A)e ‘ ^ m(8,t)dsdt (4) 

^qo «»oo 

where A is the sampling interval. 
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Putting equation 3 Into equation 2« and using equation 4 to 
Include the term, 

p(s,t)q(s,t) 

ni(S/t) ■ ■ -I ■ ■ ■ . (5) 

p(8«t)p{s,t) + n*n(s/A)ir(t/A) 

The effective power pattern, given by ¥(x,y) > I M^P^(x,y)j 
Is the Fourier transform of 

p(s,t)p(s,t)q(s,t) 

tj;(8,t) ■ ■■■■■ -■■■ - . -.I- , (6) 

p(s,t)p(s,t) + n*n(s/A)II(t/A) 

In the case Q(x,y) « 6(x)6(y), so that q(s,t) ■ 1, this 
solution is that of, for example, Schell (1969, p 573). If 
the antenna is viewed as a spatial filter, this case becomes 
what Davenport and Root (1958) call an ”opv.lmum linear filter." 
Equation 6 is a generalization of a Wiener filter (Wiener I960) . 
The role of Q(x,y) is now clearly seen: Q(x,y) is convolved 
with the Wiener filter and broadens it. 

The question that has been postponed until now is, 

A 

What is the meaning of B, the estimated brightness temperature? 
What is meant by the notion of the true brightness distribution 
must now be examined. 

One notion of the true brightness temperature 

B(x,y) might be that of the brightness temperature at the point 
(x,y) . This definition will not do, however, since 
tlTi.s can only be observed with an antenna with an infinite 
ap'^irture. Another notion is that the true value of B(x,y) 
is given by a uniformly weighted average over a finite area. 
This is equivalent to choosing Q(x,y) « 1 within that area. 
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and ■ 0 outaida tha araa. Tha apartura illumination function 
p(a,t) goat to xaro at tha adgaa of tha apartura i 
tharafora tha bahavlor of q(a»t) at larga spatial fraquancias 
cannot affact 4»(s,t) and are lest in Tharafora 

tha notion of tha trua valua of B(x«y) being a uniformly 
weighted average over a finite araa A is also not useful. 

The only useful definition of tha trua valua of B(x,y) 
is a non-uni formly weighted average over an infinite area 
(where the weights become vanishingly small suitably far 
from the direction the antenna is pointing) ; the weight 
is given exactly by H* (x/;') . The point here is that only the 
last notion of the true brightness distribution is observable 
with a finite-aperture antenna; therefore the quality of the modi- 
fied antenna pattern 4'Cx,y) must be measured. against some reali- 
zible pattern. For example, if the error is mcanured as being 
the difference between V(x,y) and a delta-function, the best 
fit will be found to be one with a neurrew main lobe, and hugh 
side lobes that alternate sign. This may be a good solution 
mathematically, but is very far from being a useful correction 
to the data. 
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b) Moia« amplification 

Th« noise amplification is found from 


a 


2 



j ) 
00 00 


{m(s»t) }^dsdt. 


(7) 


The noise amplification will be large whenever 4»(s«t)/p(s,t) 
is large; this places one constraint on what patterns can be 
synth<=jsized. One immediate consequence is that any pattern V(x#y) 
that has non-zero spatial components where p(s,t) • 0 (i*e.« 
tp(s,t)/p(8#t) «) will also have infinite noise amplification. 

The relative magnitudes of the spatial components of P(x,y) can 
be modified wherever those components are non-zero. A response 
cannot be created in those spatial components of ¥(x,y) that are 
zero in P(x,y). 

The noi«'!% amplification can be calculated for the case where 
the aperture is uniformly illuminated— p(s) « (1 - |si)H(iss) — 
and the Wiener filter form for V>(a) ia used: 


1 - |s| 

m(s) « - - — n(>n8) 

1 - Isl + n* 


Then 

1-1 1 

a* « tan '‘■(l/n) - (9) 

n 1 + n* 


(All of the examples given '^e for a one-dimensional aperture; 

a second dimension adds considerable complication but no new 

understanding.) If the Wiener filter is used, which is 

equivalent to choosing 0(x) to be a {-function, the noise amplification 

can become large— in fact, as n^-^0, although remains 

finite— and it is partly to reduce the noise amplification 

that other forms for Q(x) are sought. 
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c) A set of roalizable pattaras 

The Fourier transform relations expressed above can be 

used to investigate the properties of realizable effective 
antenna patterns. He shall develop a ccxnplete set of possible 
antenna patterns— complete in the sense that any realizable 
pattern is a linear combination of these basis patterns— r 
for simplicity, only one-dimensional apertures are considered. 

The aperture has a width of one; the autocorrelation function 
of the electric field distribution in the aperture- p(s)— 
will be zero outside the interval (-1,1). Let 'f(x) be 
the effective antenna pattern and t|i(s) its Fourier transform. 

To start with, we must realize that it is possible to 
write q(s) as 

q(s) (p(s)^ + n*) 

q'(s) » (10) 

P(s)-^ 

for some q(s),* then tt*(s) “ q(s) and m(s) ■ \)>(s) / p(s) 

* q(s) / p(s). The question then is, Hhat are the possible 
forms that q(s) might have? 

The first requirement is that m(s) must not be too large, 
since the noise eunplification is / m(s)^ds. The second 
requirement is that it be possible to write 

m(s) • &Q * I a^cos(nirAs) (11) 

where A/2 is the sampling intvrval; this requirement comes about 
from the definition of m(s) . In equation (11) , attention is 
limited to symmetric antenna patterns, so the sine te.nns do 
not appear. 

Consider, as an example, the case where A ■ 1 and 
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the aperture is uniformly illuminated, so that p(s) ■ 1 - ls|. 

Then for m(s) given in equation 11, the only aperture distributions that 
can be generated are linear combinations of 

2 (1 - |sl) • cos(nirs)II(s/2) . (12) 

The resultcuit antenna patterns are the Fourier transforms of 
<j». given by 

2 n* + 4x^ _ 

^ (x) , {1 - (-1)" cos2ttx}; (13) 

(n^- 4x^)2 

H'q (x) a sine (x) ,which is the power pattern of a uniformly 
illuminated aperture. The functions are shown in Figure 

2 for n=0 to 5. They have the property that they have maxima 

at X s n/2, and have zeroes at integral values of x (n even) 

or half-integral values of x (n odd) . Individually, the 
'C^(x) are not useful antenna patterns; the aim is to find a 
linear combination that is useful (i. e., goes to zero quickly 
for X > 1) . 

For Q(s) we now have to choose a desirable function, but 
it must be bandwidth limited. Several families of functions 
are possible; the one that we shall use as an example is 

j„(2irx) 

F (x) » 2"^Sl — (14) 

(2ttx)" 

^ V 

where jj^(x) is the n spherical Bessel function of the first kind. 

The Fourier transform of F^(x) is (1 - s^)”lI(s/2). Some other 

n 

possible families are given in Table 5. The first five functions 

F (x) are shown in Figure 3. 
n 
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The noise amplification for each pattern F^(x) (n > 0) is 


n 2^"{(2n-l) l}^ 


o* 


2n-l 


-1 


ILZJSI 

(1 - Isl) 


2n-l 

2n 


(4n -1) 1 


ds* 


(15) 


For n ' 0, the noise amplification is infinite unless t|^(s) ■ 

II(s/2b) for some b < 1. Then o*q ■ 2b/(l-b). The interesting 
feature of this family of functions is that is a decreasing 
function of n, and also the sidelobes (defined as anywhere |x| > 1) 
decrease with increasing n, although the width of the main lobe 
increases with increasing values of n. 

In order that be realizable, it must be possible to 

expand (1 - s*)”/(l - |s|) in a Fourier cosine series; i. e., 


(1 - s*)" 
(1 - |sl) 


N 

a- I a cos (nits). 
0 “ n 

n«l 


(16) 


When m(s) is specified by such a series, the noise eunplification 
can bs found by 

a* - 2a J + y a* . (17) 

The Fourier cosine series expansion coefficients for F 2 , 
and F^ are given in T£ible 6. Most of the noise eunplifica~ 
tion is due to the first term in each series. 

Since, in practice, only a finite nyunber of antenna temperatures 
can be included in the calculation, the series will have to 
be truncated after a few terms; this will not significantly 
reduce the noise eunplif ication, but will change the shape 
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of the pattern. Table 6 shows that the relative 
Importance of the higher-order terms increases as n increases; 
the price one pays for the low noise and low sidelobe level of 
the higher-order functions is that more and more 

antenna temperatures must be included in the calculation of 

and the width of the main beam increases. 

In the earlier section dealing with the matrix formulation 
of the problem, Q(x,y) was chosen— somewhat arbitrarily— 
to be a Gaussian. The main reasons for this choice were that 
a Gaussian form leads to workable formulae, and that it fits 
the author's idea of what is a "desirable” antenna pattern. 

The role of Q(x,y) is really to supply a goal toward which the 
antenna pattern is modified. When the formalism is changed to 
consider Fourier transforms, the use of a Gaussian becomes 
akward, since the Fourier transform of a truncated Gaussian 
is not readily evaluated. For this reason — mathematical 
convenience— a different form was taken for Q; i.e., 

Q * calcuations discussed here are only for 

the purpose of elucidation of a method— a method that will 
not be practical in any event for real scanning antennas— 
the exact form of Q is not im.portant. What is important 
is the realization that the anteima pattern can be modified 
by forming linear combinations of the antenna temperatures, 
and that the limitations as to what is possible — the aperture 
distributions ^ must be linear combinations of the functions 
i|»jj(s) defined in equation 12— can be clearly stated. 

Another feature that warrants comment is that explicit 

-2 

mention of the signal- to-noise ratio n disappeared after 
equation 8 to be replaced with discussion of the noise amplification 
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a*. This emphaalsQS somewhat better the Idea that* in general, 
that the modification of the antenna pattern is performed only 
with a concomitant change in the level of the noise. One must 
be careful not to adopt a method of antenna pattern modification 
that produces a beautiful pattern but has infinite noise amplification. 


i. i 
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d) Relationship between resolution and noise amplification 

The functions P^(x) = 2"‘*‘^(2irx) ”"j^(2irx) increase in 

halfwidth as n increases; the full<>widths at half*>maximum (FWHM) 
are given in Table 7 for the first seven functions (x) . 

It is not hard to show that, when the antenna pattern is 
modified, if the resolution is to increase, the noise 
£unplification must increase also. 


Define the fraction of the power contained in the 
interval (-Xq,Xq) in the main beaim, where Xq < %, to be 


M 


- f 


'?(x)dx 


i|<(s) 


-X, 


sin 2 itxqS 

2itXqS 


ds 


(18) 


(for a symmetric pattern) . The noise amplification is 


a 


2 


1 

* 

JiM 

J [ p(s) 
1 


2 

ds. 


(19) 


If h(s) is any ftmction defined on the interval (-1,1) that 
has the properties that h(0) - h(l) » 0 and h(x) > 0, 
we can replace t|;(s) with i|»(s) eh(s) , where e is a small 
positive quantity. Then, since sin(irs)/(irs) > 0 for Isj < h> 
du/de and da^/de are both > 0. This shows that the main 
beam can be narrowed only at the expense of increasing the 
noise amplification. 

There does not seem to be any simple relationship between 
the sidelobe level and the noise amplification. 
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When it comes to making a choice for Q(x,y ) , the 
desired antenna pattern# it is important to select a function 
that is not xinrealizable with the given antenna pattern. There 
is a trade-off between resolution and noise amplification; if 
Q(x#y) is chosen to be too narrow# two results will follow. 

The first is that the noise amplification will be unacceptably 
high# regardless of the sidelobe level. The second is that# 
since the fitting procedure for determining ¥(x#y) from Q(x#y) 
and the original pattern P(x#y) makes a compromise between 
making 4'(x#y) close to Q(x#y) and amplifying the noise, 
if Q(x#y) is too narrow (or otherwise ill-suited) # 
the compromise will be to have a large noise-amplification 
factor and*— often— side lobes that oscillate wildly. In 
other words# what one may get is a pattern that is not a good 
fit to Q(x#y) because Q(x) has the wrong shape. 
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e) Discussion 

The discussion that has been given here is basically 
similar to that of Schell (1969) » but the problem treated is 
significantly different; the differences have a considerable 
influence on the nature of the solution. 

The first difference is that we are concerned here with the 
problem that land-ocean boundaries— with a contrast of the 
order of 100 K— exist and the accuracy of the measurements of bright- 
ness temperature must, he considerably better than IK. If the effects 

of the antenna side lobes are not reduced to a vary low level# ob- 
servations over the ocean near land will be useless. Therefore# the goal 
here is not to improve the spatial resolution in the main lobe# 
but to reduce the level of the side lobes to a minimum. 

A second difference is that the data rate from instruments 
like SMMR and LAMMR is very high; there is not time enough 
to Fourier transform the data in real time. At the outset# 
the limitation was adopted that the estimated brightness 
temperature must be a linear combination of the antenna 
temperatures with coefficients that are independent of the 
antenna temperatures. This leads to the notion of choosing 
the coefficients so as to oroduce a certain emtenna pattern, 
rather than to minimize the error in the brightness temperatures. 

A third difference is that# while Schell (1969) seems to 
be interested in cases where the signal-to-noise ratio is 
high# in our case the signal-to-noise ratio is not so high 
that significant amplification of the noise can be tolerated. 

Therefore care has been taken to evaluate the noise amplification 
factors; some widening of the main beam may be tolerated if 
it leads to an increase in the signal-to-noise ratio. 
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£) Pattern matching at higher frequencies 

The problem that arises when the same antenna la used 
to make observations at more than one frequency Is that, at 
the higher frequency, the antenna pattern must be modified so 
that It matches the pattern at the lowest freq<jency. Let and 
be two frequencies, with < Vj. First, the effective 
antenna pattern— Is found at the lowest frequency. Then, for 
the higher frequency. It Is necessary to find a set of coef- 
flcents so that the effective pattern at the higher frequency 


ij/jCs) » t|/j^(s) 


( 20 ) 


or 

mjCs) « iJ/j^(s)/p2(8) . (21) 

The coefficients for the higher frequency are found from the 
Fourier transform of m2(s). Since P2(s) will not be zero 
anywhere where i(>^(s) Is non-zero, this ratio will not be large 
2mywhere . 


k: 
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V. CONCLUSIONS 

The calculations reported In this paper Illustrate a method 
for processing the data from microwave remote-sensing instruments 
so as to modify the antenna pattern. The main objective in 
modifying the antenna pattern is to suppress the response 
in the far sidelobes. The pattern chosen for study is not 
entirely realistic in that it i) does not have distinct 
side lobes, and ii) it is azimuthally symmetric. The pattern 
was chosen so that the values of and could be expres- 
sed analytically. 

Given a real antenna pattern, these coefficients would 
have to be calculated by numerical integration. Finding 
the patterns for various combinations of a^ , rf^ , and Ax 
would show how that particular pattern could be modified 
to adveuitage. While the details would differ from the 
results given here, the following general features should 
still be true: 

1) It is possible to smooth the antenna temperatures in 

such a way that the noise is reduced, but only 
the main lobe is broadened, while the level of 
the side lobes is reduced. 

2) Up to a point, it is possible the reduce the level 

of the side lobes without broadening the main lobe; 
this can be done with a reasonable loss due to 
amplification of the receiver noise. 


Nature puts a limit on the resolving power of a finite 
aperture; success in modifying the pattern of a single-dish 
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antenna depends on not asking the impossible. It is not pos- 
sible to get an estimate of the brightness at a point that 
is unaffected by the brightness distribution at surrounding 
points; if we give in to the inevitable and ask for an estimate 
of the brightness averaged over an area that is not too small, 
the quality of the results is improved significantly. This 
can be seen by comparing Figures 1-e and 1-f; the only difference 
is in the width of the Gaussian used for Q. The price of trying 
to increase the resolution is a drastically increased level 
in the side lobes and a huge increase in the no:>'3e amplification. 
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TABLE CAPTIONS 

Table 1. The nolle amplification (a‘) and lidelobe levels 
for different effective antenna patterns. The widths 
are given for the -10 db and -20 db sidelobes; the level 
of the sidelobes are given where there are separate 
sidelobes. For all of these calculations, Ax « 1. 

Table 2. The same information as Table 1, for Ax ■ 2. 

Table 3. The full width at half-maximum (FHHM) for each 
beam and the distance at which the contribution to 

the antenna temperature falls off by a factor of 1000. For 
all calculations, Ax > 1. 

Table 4. The same information as Table 3, but Ax ■ 2. 

Table 5. Some bandwidth- limited fvmctions and their Fourier 
transforms. 

Table 6. Series expansions of 

Table 7. Hals-widths fo the functions F. (x). 

n 
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Ax - 1 



S/N 

a* 

1.2 

10 

0.068 

1.2 

100 

0.18 

1.2 

1000 

0.32 

1.2 

OS 

0.56 

1.0 

10 

0.079 

1.0 

100 

0.24 

1.0 

1000 

0.50 

1.0 

00 

1.58 

0.8 

10 

0.093 

0.8 

100 

0.32 

0.8 

1000 

0.97 

0.8 

00 

7.86 

0.5 

10 

0.12 

0.5 

100 

0.55 

0.5 

1000 

2.29 

0.5 

00 

235. 


TABLE 1 


X ■ 

-10 db 

X ■ 

-20 db 

sidalobe 

levol 

3.4 

4.7 



2.7 

3.4 

21.5 


2.4 

3.2 

25.0 


2.4 

3.3 



3.3 

4.4 



2.6 

3.1, 4.1 

19.2 


CM 

• 

CM 

2.8 

24.3, 

24.0 

2.1 

3.1 



3.2 

4.2 



2.4 

2.9, 4.1 

17.5, 

25.7 

2.1 

2.5, 3.3 

18.9, 

22.3 

1.9 

2.8 

28.7, 

29.2 

2.9 

3.9 

25.7 


2.2 

2.6, 4.0 

15.4, 

22.4 

1.8 

2.1, 3.2, 
4.2 

14.9, 

29.3 

19.6 

1.5 

3.0, 4.0 

18.2, 

20.1 

19.5 


TABLE 2 
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Ax B 2 
a2 

S/N 


X ■ 

-10 db 

X ■ 

-20 


side lobe 
level 


1.2 

10 

0.29 

3.4 

4.3 


26.7 



1.2 

100 

HI 

• 

o 

2.8 

3.6, 

5.1 

17.6, 

23.2 


1.2 

1000 

0.99 

2.5 

3.5, 

4.6 

19.8, 

24.3 


1.2 

00 

1.05 

2.4 

3.4 


21.1, 

25.3 


1.0 

10 

0.34 

3.3 

4.2 


23.6 



1.0 

X‘; * 

0.94 

2.6 

3.5, 

5.0 

18.3, 

24.3 


1.0 

1000 

1.39 

2.2 

3.1 


24.7, 

29.4 


1.0 

00 

1.50 

2.1 

3.0 





0.8 

10 

0.41 

3.2 

4.1 


22.2 



0.8 

100 

1.29 

2.3 

3.3 


20.0, 

26.3 


0.8 

1000 

2.06 

2.0 

2.5 


20.2, 

22.3, 

26.5 

o 

• 

00 

00 

2.25 

1.9 

2.3, 

3.2 

18.1, 

19.9, 

24.2 


Original Beam 

2.4 

4.1 

Gaussian a^ = 1 

2.1 

3.1 
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TABLE 3 


Ax » 1 




a * 

S/N 

FWHM 

^000 

1.2 

10 

3.8 

2.9 

1.2 

100 

3.6 

2.4 

1.2 

1000 

2.8 

2.2 

1.2 

00 

2.6 

1.7 

1.0 

10 

3.7 

2.5 

1.0 

100 

3.0 

2.7 

1.0 

1000 

2.6 

2.1 

1.0 

00 

2.4 

1.5 

0.8 

10 

3.5 

2.0 

o 

• 

00 

100 

2.9 

2.8 

oo 

• 

o 

1000 

2.5 

2.5 

0.8 

00 

2.2 

1.4 

0.5 

10 

3.4 

2.2 

0.5 

100 

2.7 

2.9 

in 

« 

o 

1000 

2.0 

3.2 

0.5 

00 

1.7 

3.0 



TABLE 4 


Ax > 2 


2 

S/N 

FWHM 

^000 

1.2 

10 

3.4 

2.6 

1.2 

100 

2.7 

2.9 

1.2 

1000 

2.5 

2.2 

1.2 

« 

m 

• 

CM 

2.2 

1.0 

10 

3.2 

2.5 

1.0 

100 

2.5 

2.3 

1.0 

1000 

2.3 

1.5 

1.0 

00 

2.3 

1.5 

00 

• 

o 

10 

3.0 

2.6 

0.8 

100 

2.6 

2.1 

o 

• 

00 

1000 

2.3 

2.2 

0.8 

00 

2.2 

GO 

• 

CM 

Original Beam 

2.4 

5.9 


Table 5. some bandwidth-limited functions and their Fourier transforms 


Function 
(1 - s*)”n(%s) 

(1 - s*)"^‘*n(Jjs) 

(-i)"2T^(s)n(i|s) 
(I - 

(-i)"p^(8)n(%s) 
ccs^'cJjirs) n (Hs) 


Transform Remarks 

j (x) s spherical Bessel function 
of the first kind. 

J = Bessel function of the 
first kind. 

J (2wx) T = Chebyshev polynomial of 

" " the first kind. 

^ ~ Legendre polynomial. 

2^""r(ns-l) 


n + l J_(2lfX) 
- 2 - 


(2irx) 


n 


^ ! — 

(2iix)" 


r(l+ijn+2x)r(l+»in-2x) 


Table 

6. Series 

expansion 

of F = 2”^^ n! j„ 
n n 

(x)/x”. The 

values 

of a_ are 
n 

given 

separately 

for n»0, n 

1 even, and n odd. 

^ n • 

a_ cos (nss) . 
n 


^1 


^2 

^3 



n = 0 

3/2 


11/12 

7/10 




A 


-10/(nTT)^ 

- 2 

168 


SI 


(nw) ^ 

(nir) ^ 


n odd 

- 4 


6 24 

- 2 

120 

480 

(nir) ^ 


(nir)^ (nir)^ 

2 * 
(nn)^ 

(nn) ^ 

(nr)® 


a 


2 


n 


14/3 


198/105 


1486/llb5 
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Table 7. Half-widths of the functions F (x) 

n 


2»+lnl 


vr* 

<ltx)“ 


n FWHM 

0 0.60 

1 0.80 

2 0.95 

3 1.09 

4 1.22 

5 1.32 

6 1.42 


FIGURE CAPTIONS 


Figure 1. The original emtenna pattern and a Gaussian of 
unit width (Fig 1-a) . Computed effective patterns 
are shown in 1-b to 1-g. The patterns are normalized so that 
the integral is unity. 

Figure 2. The functions n = 0, . . 5. 

n+1 3n^2irx) 

Figxire 3. The functions 2" n! — - t n = 0, « . 4. 

(2irx) 
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Figure 1-a. 
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Figure 1-d. 



Figure 1-g. 
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Figure 1-b 


Figure 1-c 
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